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Abstract 

We demonstrate strong coupling between the flexural vibration modes of a clamped-clamped mi- 
cromechanical resonator vibrating at low amplitudes. This coupling enables the direct measurement 
of the frequency response via amplitude- and phase modulation schemes using the fundamental 
mode as a mechanical detector. In the linear regime, a frequency shift of 0.8 Hz is observed for a 
mode with a line width of 5.8 Hz in vacuum. The measured response is well-described by the ana- 
lytical model based on the Euler-Bernoulli beam including tension. Calculations predict an upper 
limit for the room-temperature Q-factor of 4.5 x 10 5 for our top-down fabricated micromechanical 
beam resonators. 
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Nonlinear interactions between the vibration modes in micro- and nanomechanical res- 
onators have attracted significant interest recently. In extensional structures, such as 
clamped-clamped bridges, the modes are coupled by the displacement-induced tension 
which yields a quadratic relation between the resonance frequency of the mode considered 
and the amplitudes of the other modes. In singly-clamped cantilevers the displacement- 
induced tension is absent; here the inextensionality condition couples the horizontal and 
vertical displacements of all modes, resulting in qualitatively similar dynamics {2 ^J. Sev- 
eral applications and consequences have been put forward based on the modal interactions, 
such as enhancement of the dynamic range modification of the resonator damping by 
employing phonon-phonon cavities , frequency stabilization |8j, the study of relax- 
ation mechanisms {9] and linear frequency conversion [lO]. In thermal equilibrium, via the 
modal interactions the displacement fluctuations in one mode give rise to frequency fluc- 
tuations in the other modes, thus broadening their spectral line. In recent theoretic work 
these frequency fluctuations were quantified for a carbon nanotube [llj], yielding a bound- 
ary for the experimental Q-factor of such resonators at room temperature. Experiments 
on suspended carbon nanotube resonators in the Coulomb blockade regime demonstrate 
that single-electron-tunneling processes provide a strong electrostatic coupling between the 



modes, in addition to the mechanical mode coupling [12|, 

While the recent experimental work has focused on the mode coupling in strongly 
driven resonators exhibiting nonlinear vibrations, in this Letter we investigate these inter- 
actions in the low-amplitude regime. We demonstrate that the modal interactions play a 
significant role in the dynamic behavior of a linear resonator, as the vibrations of a weakly 
driven mode modulate the motion of a second vibration mode. We employ this coupling 
to perform swept-frequency type measurements of the linear frequency response of a high 
frequency (target) mode, by measuring the induced amplitude (AM) and phase modulation 
(PM) in a low frequency (probe) mode which is weakly driven at a fixed frequency. This 
provides a practical way to measure the frequency response, and it should be contrasted 
to the scheme presented earlier where the frequency response of the target mode is re- 
constructed from a series of resonance frequency measurements on the probe mode. The 
observed modulation depth is in agreement with a model based on the Euler Bernoulli beam 
including tension Using this model with the parameters extracted from experiment, 
the frequency broadening of the fundamental mode that occur when the high amplitude 
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mode is in thermal equilibrium is estimated. An upper limit Q = 4.5 • 10 5 is found for our 
micrometer-scale silicon beam resonator. This indicates that modal interactions play a sig- 
nificant role in the experimentally observed Q-factors of micromechanical resonators. 

To detect the modal interactions in the linear regime we deploy the fundamental flex- 
ural mode as the probe, by driving it at a fixed frequency f pro be- Its response is modulated 
by the averaged tension induced by the displacement of a second mode of the resonator, 
which is driven at ftarget- Fig- 1(a) shows the principle of this AM scheme. The experiments 
are conducted using silicon beams with dimensions L x w x h = 1000 x 2 x 2 /im 3 , fabri- 
cated by anisotropic wet etching of silicon-on-insulator wafers with an (llO)-oriented device 

layer 14|. Fig. 1(b) shows the resonator and the setup; the beam is driven by applying an 

I I 

alternating current in the presence of a permanent magnetic field [15| , via a thin conductive 
layer evaporated on top (70 nm of chromium/gold). The driving force and the detector 
couple to the in-plane motion of the device. To eliminate broadening of the resonance peak 
by the viscous force, the experiment is conducted at a pressure of ~ 10 -4 mbar. Frequency 
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FIG. 1: (a) Amplitude modulation scheme to probe the modal interactions at the low amplitudes. 
The response of the probe mode is modulated by the displacement-induced tension of the target 
mode, (b) Measurement setup and colored scanning electron micrograph (top view) of the silicon 
(110) beam resonator. The motion is in the plane of the figure, (c) Linear frequency response of 
the probe and target modes. 
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response measurements when one mode is driven and the other mode is in thermal equi- 
librium are shown in Fig. 1(c). The linear resonance frequencies for the first and third 
mode (n = 1,3), f^i = 26276.6Hz and fR j3 = 127294Hz, and the corresponding Q-factors, 
Qi = 4527 and Q3 = 7476, are obtained from damped-driven harmonic oscillator fits (solid 
lines) l^l . 




The sensitivity of the probe mode to the motion of the target mode varies with f pro be 
and maximizes when the second derivative of the probe mode frequency response is zero. In 
Fig. 2(a) the fitted amplitude (left) and phase (right) responses of the fundamental mode 
(grey) of Figure 1 are repeated, together with their second derivatives (blue solid lines). 
The inflection points, marked B and D, occur at detuning Af m ^7=, whereas the phase 
modulation maximizes on resonance, Af = (C). Figure 2(b) shows measurements of the 
magnitude (left) and phase (right) response of the probe mode while sweeping f tar get through 
the resonance frequency of mode 3. The (fixed) drive frequencies correspond to A-E in panel 
(a). At fprobe = 26.76kHz, (A), the probe mode is off-resonance and the target mode is not 
detected. In B, on the positive slope of the frequency response curve of the probe mode, the 
amplitude modulation maximizes. The tension of the target mode tunes the resonance fre- 
quency of the probe mode to a higher value, and produces in a dip in the amplitude response 
of the probe mode. When the probe mode is on resonance (C) the amplitude modulation 
is close to zero: here the phase response maximizes. Probing at a higher frequency results 
in maximization of the AM response in D, and off-resonance (E) the modulation is zero. 
These measurements demonstrate that the fundamental mode can be conveniently used to 
probe the response of a mode that is 5 times higher in frequency, in a swept-frequency type 
measurement. 

Figure 2(c) shows the full measurement from which the traces are extracted. The 
target mode is swept along the horizontal (fast) axis, while the frequency of the probe mode 
is incremented along the vertical (slow) axis. The color scales represent the amplitude (left) 
and phase (right) response of the probe mode. The probe response when the target mode is 
off- resonance is subtracted to reveal the resonant features more clearly. Figure 2(d) shows 
the modulation depth as a function of the probe mode frequency, when the target mode is 
driven on-resonance. The AM modulation depth is 10 mV, which corresponds to a frequency 
tuning of 0.8 Hz. This frequency shift is of the same order as the line-width of the resonator, 
5.8 Hz, indicating the strong coupling between the modes vibrating in the linear regime. 
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Notably, in these measurements the applied drive force for both modes is the same as in 
Fig. 1(c), and no apparent sign of nonlinearity is observed in the line shapes. In Fig. 2 
(d) the back-action from the target mode causes a small but noticeable upwards shift of the 
resonance frequency of the probe mode response. 

To describe the modulation between the two linear driven modes, we calculate the 
complex response of mode 1 while driving mode 3 through resonance, using the model de- 
veloped in Ref. [lj]. The amplitudes ai j3 of modes 1,3 driven at frequency are then given 
by: 

( w fU - w i + iuiU R>1 /Q 1 + Jail 2 /!! + ^{\a 3 \ 2 I u h3 + N^is)) ! = J /ac&dx, (1) 

( w h,3 ~ u 3 + iu&>Rfl/Qs + ^|a 3 | 2 / 2 3 + ^(|a 3 | 2 /n/ 33 + |a 3 | 2 / 2 3 ))ai = J f^&dx. (2) 

The scaled Lorentz force per unit length acting on the beam is f aCji = j^p -Bii, where 
I\ = 715 pA and 7 3 = 10.7 nA are the rms currents through the resonator and B = 1.9 T 
is the permanent magnetic field. For a stress-free beam, the integrals Iy = J ^(a;)^'-(x)dx 
are In = 12.3, J 33 = 98.9 and ii 3 = —9.7, and the average displacements per unit deflection 
are Jq 1 £idx' = 0.83 and Jq 1 ^dx' = 0.36 for the beam-like mode shapes Young's modulus 
equals E = 169 GPa for our device, which bends about the silicon (110) axis Figure 2(d) 
(solid lines) shows the calculated modulation depth as calculated by numerically solving the 
coupled equations of motion. Agreement is found between the measured and the calculated 
sensitivity curves, where the phase shift is quantitatively predicted by the model, and the 
magnitude is scaled by a constant factor as the detector gain is not accurately known in this 
experiment. 

The above analysis demonstrates that the modal interactions give rise to a significant 
tuning of the modes in the low-amplitude regime. The driven displacements of mode 1 in 
the experiment exceed the thermomechanical fluctuations by only 2 orders of magnitude, 
and introduce a significant frequency tuning when compared to the resonator bandwidth. 
It is interesting to consider the frequency shifts when the modes are in thermal equilib- 
rium and the external driving force is zero. Here, the displacement fluctuations give rise 
to frequency fluctuations via the displacement-induced tension. For carbon nanotube res- 
onators this broadening was calculated using a thermodynamic approach and a discretized 



toy model of the resonator 11[. With the continuous Euler-Bernoulli theory, which accu- 



rately describes the experimentally observed interaction between the modes at low ampli- 



tudes, the frequency shift for our top-down fabricated micromechanical resonators equals 



0.09Hz/nm 2 . For an rms thermal displacement a3 jt h = y keT/mo;^ 3 = 0.026 nm this figure 
implies an upper bound for the experimentally observed Q-factor of mode 1 of approximately 
Qi,max < 4.5 • 10 5 at room temperature, due to the interaction with mode 3. This value will 
be lower when the motion of other degrees of freedom are taken into account. 

The predicted upper bound is lower than the experimental Q-factors for string 
resonators under tension, for which values over a million have been reported at room- 
temperature [18j. The difference can be explained as follows. The string devices incorporate 
a large residual tension, which forms the main contribution to the restoring force 19[. As a 
result, the flexural displacements due to thermal energy are small, and so is the displacement- 
induced tension. The dispersion due to a fluctuating displacement via the modal interactions 
therefore reduces with the residual tension. This is mathematically expressed via the in- 
teraction matrix ly, which for devices with negligible bending rigidity (strings) contains 
only nonzero diagonal terms (I^j = 0), thus limiting the available dispersive interactions 
for strings. Finally, in a perfect string the spectrum is harmonic with f^ n = nfj^i. In this 
case the energy is recycled between the modes as the motion of one mode gives rise to a 
parametric excitation of the other modes via the displacement-induced tension. In a perfect 
string this parametric excitation is effective, as it always occurs at integer multiples of the 
resonance frequencies and all harmonics are parametrically degenerate. 

In conclusion, we have studied the modal interactions in clamped-clamped resonators 
in the regime of linear vibrations. At drive strengths that are 5 orders of magnitude lower 

n 

than in previous experiments on micromechanical resonators [1] , we observe a strong mechan- 
ical interaction between the modes, inducing frequency shifts on the order of the resonance 
line width in vacuum. The interaction is employed to perform swept-frequency type mea- 
surements employing by measuring amplitude and phase modulation of the fundamental by 
a mode that is 5 times higher in frequency. The presented scheme, i.e. coupling the motion 
to a mechanical detector, may be favorable compared to an electronic detector as it does 
not suffer from signal loss due to parasitic capacitance. The amplitude and phase modula- 
tion is quantitatively described by an Euler-Bernoulli beam model including tension. This 
model predicts that the fluctuating displacements in thermal equilibrium result in frequency 
broadening, and limit the room-temperature Q-factor of our micromechanical resonators to 
4.5 • 10 5 . In reality this value will be lower since the model only considers broadening by the 
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third mode, and neglects the contribution of the other modes. 
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FIG. 2: (a) Frequency response of the probe (fundamental) mode and its second derivative, (b) 
Amplitude (left) and phase (right) response of the probe mode driven at the frequencies indicated 
in (a), while driving the target mode through resonance. (c) Amplitude (color scale, left) and phase 
(right) of the probe mode, after subtracting the harmonic oscillator response background, (d) 
Sensitivity of the probe mode in AM and PM modulation schemes as a function of the drive 
frequency of the probe (solid lines: calculation). The target mode is driven on resonance, (e) 
Response as calculated by solving the coupled equations of motion (see text). 
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